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Abstract. We study the dynamics of a certain discrete model of interacting particles that 
comes from the so called shuffling algorithm for sampling a random tiling of an Aztec diamond. 
It turns out that the transition probabilities have a particularly convenient determinantal 
form. An analogous formula in a continuous setting has recently been obtained by Jon Warren 
studying certain model of interlacing Brownian motions which can be used to construct Dyson's 
non-intersecting Brownian motion. 

We conjecture that Warren's model can be recovered as a scaling limit of our discrete model 
and prove some partial results in this direction. As an application to one of these results we 
use it to rederive the known result that random tilings of an Aztec diamond, suitably rescaled 
near a turning point, converge to the GUE minor process. 



1. Introduction 

There has been a lot of work in recent years connecting tilings of various planar regions with 
random matrices. One particular model that has been intensely studied is domino tilings of a 
so called Aztec diamond. One way to analysing that model, [IHS], is to define a particle process 
corresponding to the tilings so that uniform measure on all tilings induces some measure on this 
particle process. 

In this article we will study the so called shuffling algorithm, described in [51lllj. which in 
various variants can be used either to count or to enumerate all tilings of the Aztec diamond or 
to sample a random such tiling. 

The sampling of a random tiling by this method is an iterative process. Starting with a tiling 
of an order n — 1 Aztec diamond, a certain procedure is performed, producing a random tiling 
of order n. This procedure is usually described in terms of the dominoes which should be moved 
and created according to a certain procedure. We will instead look at this algorithm as a certain 
dynamics on the particle process mentioned above. 

The detailed dynamics of the particle process will be presented in section [2] and how it is 
obtained from the traditional formulation of the shuffling algorithm is presented in section [U 
For now, consider a process x{t) — (X 1 (t) , . . . , X m (t)) for t = 0, 1, 2, . . . , where X k (t) = 
(Xi(t), . . . , X%(t)) £ Z k . The quantity X\ (t) represents the position of the i:th particle on line 
j after t — j steps of the shuffling algorithm have been performed. (The reason for the t — j is 
technical convenience.) We will show that 

Theorem 1.1. For fixed k, consider only the component X k {t) from x{t) rescaled according to 

(1) x?{t) 



X?{Nt) -\Nt 



and defined by linear interpolation for non-integer values of Nt. The process X n (t) converges 
to a Dyson Brownian motion with all particles started at the origin as N — > oo, in the sense of 
convergence of finite dimensional distributions. 

The full process (-t(t))t=o,i,... nas remarkable similarities to, and is we believe a discretization 
of, a process studied recently by Warren, [12] . It consists of many interlaced Dyson Brownian 
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motions and is here briefly described in section[3] We will denote that process (X(t)) t >o. There 
is reason to believe the following. 



Conjecture. Consider the process (x (i))t=o,i,... rescaled according to 

(2) xm 



and defined by linear interpolation for non-integer values of Nt. The process x(f) converges to 
Warren's process X(t) as N — > oo, in the sense of convergence of finite dimensional distributions. 

The key to our asymptotic analysis of the shuffling algorithm is that the transition probabilities 
of (X k , X k+1 ) can be written down in a convenient determinantal form, see proposition ^. 21 These 
formulas mirror beautifully formulas obtained by Warren. 

As an application of our results we will use it to rederive an asymptotic result about random 
tilings near the point where the arctic circle touches the edge of the diamond. This result was 
first stated in [6] and proved in [4j. 

Recall that the Gaussian Unitary Ensemble, or GUE for short, is a probability measure on 
Hermitian matrices with density Z~ 1 e~ TtH I" 1 where Z n is a normalisation constant that depends 
on the dimension n of the matrix. Let H — (/i rs )i<r,s<n a GUE matrix and denote its principal 
minors by Hj = (h rs )i< r ,s<j- Let A J = (Aj, . . . , Aj) be the eigenvalues of Hj. Then A = 
(A 1 , . . . , A") e M«(«+i)/2 i s the so called GUE minor process. 

Theorem 1.2 (Theorem 1.5 in [4 .). Let the l"( n+1 )' ' 2 -valued process x(t) = (X 1 ^), . . .,X n (t)) 
be a rescaled version ofx(t) with 



o) xm = —t ] 
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Then x (t) — > A as t — > oo in the sense of weak convergence of probability measures. 

To put this in perspective, let us note that a similar result for lozenge tilings is known from 
Okounkov and Reshetikhin [TU]. They discuss the fact that, for quite general regions, that close 
to a so called turning point the GUE minor process can be obtained in a limit. A turning point 
is, just as in our situation, where the disordered region is tangent to the domain boundary. 



2. The Aztec Diamond Particle Process 

We will here content ourselves with stating the rules of the particle dynamics that we will 
study. The reader will in section |4] find a description the traditional formulation of the shuffling 
algorithm and how that relates to the formulas below. 

Consider the process (x(t)) = (A 1 ^), . . . , X n (t)) for t = 0, 1, 2, where X k {t) = 
(X k (t), . . . , X*(t)) e Z fc . It satisfies the initial condition 

(4) X k (0) = x k 

where = i for 1 < i < j. At each time t the process fulfils the interlacing condition 



(5) 



x k (t)<xt\t)<x k +1 (t) 
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and evolves in time according to 

x}(t)=xt(t-i)+pi(t) 
x{(t)=xi{t-i)+p{(t) 

- l{X((t - 1) + /3{(t) = X{-\t) + 1} for j > 2 
X>(t)=X 3 ;(t-l)+0>(t) 



(6) 



l{X{(t - 1) + (3{(t) = X]Z\{t)} for j > 2 



Xj(t)=Xj(t-l)+0>(t) 

-i{xi(t-i)+ ( 3i(t) = xr i (t) + i} 

+ l{Xi(t-l) + Pl (t) = XlZl (t)} for j > 3 and 1< i < j. 

for t = 1, 2, . . . where all the j t are i.i.d. unbiased coin tosses, satisfying P[f3\ (1) = 0] = 

P[/3 1 1 (l) = l] = i. 

One way to think about this is that at each time t, this is a set of particles on n lines. The k:th 
line has k particles on it at positions Xf, . . . , X%. At each time step each of these particles either 
stays or jumps one unit step forward independent of all others except that the particles on line 
k can force particles on line k + 1 to jump or to stay to enforce the the interlacing condition (J5]). 
Also note that the interlacing implies that Xf < X^ +1 at each time t, i.e. two particles cannot 
occupy the same space at the same time. 

As mentioned we can write down transition probabilities for this process on a particularly 
convenient determinantal form. Define Si : Z — > Z such that 5i(x) = 1 if i = x and Si(x) = 
otherwise. Let us first introduce some notation. 

(4> * ip)(x) = <p(s)ip(t) (Convolution product) 

s+t—x 

(o) ^o 

=<£("-!)* $ for n = 1,2,... 

A<fi = (Sq — Si) * (j) (Backward difference) 



V=—oo 

Acj) = (—Sq + S-i) * 4> (Forward difference) 



x-l 



A- 1 cf>(x)= <KV) 



y=-oo 



Let w n+1 > n 

w n+l,n and t = 

(7) 

where 



= {(x,y) : xi < yi < x 2 < 
0, 1, ... , define 



< 



?((x,y),(x',y , )) = det 



< x n+1 } C Z" +1 x Z". For {x,y),{x',y') £ 



At(x,x') B t (x,y'j 
C t (y,x') D t (y,y') 



• A t (x, x') is an (n + 1) x (n + l)-matrix where element (i, j) is (^'{x^ — Xj), 

• B t (x,y') is an (n + 1) x (n)-matrix where element (i, j) is A^ 1 (j)^(y' i — Xj) — l{j > i}, 

• Ct(y, x') is an n x (n + l)-matrix where element is A^*-*)^ — Xj) and 
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• D t (y,y') is an n x n-matrix where element is (^'{y'i — Uj)- 
Let W n = {x : xi < x 2 < • ■ ■ < x n } C Z n and for x G W" let 

(8) ftn(af) = JJfe - 

i<j 

Finally, after all this notation, we can state a result. 

Theorem 2.1. TTie transition probabilities of (X k ,X k+1 ) from the process x above are 

(9) !?' + ((x, y), (x 1 , y' j) := ^7T? t fe ((*, y), (*',</)) 
i/iai is 

(10) F[(X k+1 (s + t),X k (s + t)) = (x',y');(^ &+1 (s),^(^) - (x,y)] = ?t fel+ ((a;,y),(^,y'))- 

A proof is given in section E] and the reason I defined q k as opposed to defining q^' + directly 
will become obvious in the next section. 

Given the exact expressions above it is a very straightforward computation to integrate out 
the x component in expression ©. We find that the transition probabilities of (X k ) from the 
process x above is 

(11) Ft (y,y):=j^ Pt (y,y) 

where p k {y,y') ■= D t (y,y') given above. We recognise this as the transition probability for 
random walks conditioned never to intersect, a fact that is so important we state it properly. 

Corollary 2.2. The component X k (t) of x(t) is a discrete Dyson Brownian motion of k particles 
started at x k . 

This fits nicely with theorem 11.11 The component X k from x simply k simple symmetric 
random walks conditioned never to intersect, their limit is k Brownian motions conditioned 
never to intersect, which is exactly what X k from Warren's process X is. 

3. Interlacing Brownian motions 

We will now digress a bit and summarise Warren's work in [T2], so as to see the similarities 
between his continuous process and our discrete process. The reader is referred to that reference 
for more details of the construction. Consider an R" +1 x R n - valued stochastic process (Q(i))t>o = 
(X(t),Y(t)) t >o satisfying an interlacing condition 

(12) X^t) < Yi(t) < X 2 {t) < ■ ■■ < Y n (t) < X n+1 (t), 
and equations 

(13) Yi(t) =y i +f3 l {t At), 

(14) Xi(t) = y t + 7i (i A r) + Lr(t A r) - Lf {t A r) 
where 

(0i)f=i and are independent Brownian motions, 

r = inf{t > : Yiif) = Yi + i for some i}, 

r + 

J n+l 



L x = L+ +1 = and 



(15) L+(t) = / l(X t (s) = Yi(a))dL+(a) Lr(t) = / l(X t (s) = Y^s)) dLr ( s ) 
Jo Jo 

are twice the semimartingale local times at zero of Xi — Yi and Xi — Yi-\ respectively. 
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(16) 



This process can be constructed by first constructing the Brownian motions (3i and 7; and 
then using Skorokhod's construction to push Xi up from Yi_i and down from Yi. The process is 
killed when r is reached, i.e. when two of the Yi meet. 

Warren then goes on to show that the transition densities of this process have a determinantal 
form similar to what we have seen in the previous section. Let ft(x) = {2iit)^ 1 / 2 e^ x / 2 * and 
Mx) = Coo Mv) d V- Let W n < n+1 = {(x,y) e R" x R n+1 : x x < y x < x 2 < ■ • • < y» < x n+1 }. 

Define q™((x,y), (x',y')) for (x,y), (x',y') € w rhn+1 and t > to be the determinant of the 
matrix 

~A t (x,x') B t (x,y')~ 
C t {y,x') D t (y,y')_ 

where 

A t (x,x') is an (n + 1) x (n + l)-matrix where element is <fit{x'i — Xj), 
B t {x,y') is an (n + 1) x (n)-matrix where element («, j) is $t(y- — aJj) — 1(J > i), 
Ct(y,x') is an n X (n + l)-matrix where element («, j) is ^(j/^ — £j) and 
D t (y,y') is an n x n-matrix where element is <Pt{y[ — yj)- 

Proposition 3.1 (Prop 2 in [12] ) . The process (X, Y) killed at time t has transition densities 
q™, that is 

(17) <j?((x,v), {x,y))dx'dy' = V*>v[X t € dx',Y t € dx';t < r] 

Warren goes on to condition the Yj, not to intersect via so called the Doob /i-transform. The 
transition densities for the transformed process are given in terms of the those for the killed 
process by 

(18) q ?' + ((x,y),(x\y')) = ^f!- q n(x,y),(x',y , )y 

He also shows that you can start all the Xi and Yi of the transformed process at the origin by 
giving a so called entrance law, 



(19) v»(x,y) := /2 exp j-^> 2 /(2t) j jjjfo -a*) j |JJ( W , 

that is, showing (lemma 4 of [T^) that this expression satisfies 

(20) v? +s (x',y')= [ v«{x,y)q?' + {{x,y),{x',y'))dxdy. 

It is possible to integrate out the X components in that transition density and entrance law. 
The result is transition density 



(21) p?< + (y,y'):= -^{det D t (y,y') 

Kv) 

and entrance law 
(22) 



(V) ■= ^- 2 / 2 exp j-^ 2 /(2*)j Sn { y 3 ,-yA 



Now comes the interesting part. Let K be the cone of points x = (x 1 , . . . , x n ) where x k 
(xi, ...^'Jel 1 . Warren defines a process X(t) taking values in K such that 

(23) X*(t) = x* + 7 * (t) + L*'"(t) - L*' + (i) 
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where the ("f k )i, k are independent Brownian motions and L k ' + and L k ' + are continuous, increas- 
ing processes growing only when X k (t) — X^ 1 ^) and X k (t) = X k _ x (t) respectively and the 
special cases L^' + and L k '~ are identically zero for all k. 

Think of this as essentially n(n + l)/2 particles performing independent Brownian motions 
except that the k particles in X k can push the particles in X k+1 up or down to enforce the 
interlacing condition that the whole process should stay in K. 

This full process process can be constructed inductively as follows. 

(1) The process (X k ) has transition densities p k ' + and entrance law fi k . 

(2) The process (X k , X k+1 ) has transition densities q t ' + and entrance law v k . 

(3) For k = 2, . . . , n — 1 the process (X k+1 ) is conditionally independent of (X 1 , . . . , X k ~ 1 ) 
given (X k ). 

(4) This implies (by some explicit calculations) that (X k+1 ) has transition densities p k+1:+ 
and entrance law /x t +1 . 

This argument shows that the following. 

Proposition 3.2 (Warren). There exists such a process X(t) started at the origin and it satisfies 
that for k = l, . . . , n— 1, the process (X k , X k+1 ) has entrance law and transition probabilities 
q k >+. 

It is this process X that is the continuous analog of our discrete process x . 

4. Shuffling algorithm 

We will now show how relate some well known facts about sampling random tilings of an 
Aztec diamond before showing how to get the particle dynamics in section [5] 

The Aztec diamond of order n, denoted A n , is an area in the plane that is the union of those 
lattice squares [a, a + 1] x [6, b + 1] C K that are entirely contained in {\x\ + \y\ < n + 1}. A n 
can be tiled in 2™( n+1 )/ 2 ways by dominoes of size 2x1. We will be interested picking a random 
tiling. By random tiling in this article we will always mean that all possible tilings given the 
same probability. 

A key ingredient of almost all results concerning tilings of this shape is the realization that 
one can distinguish four kinds of dominoes present in a typical tiling. The obvious distinction 
to the casual observer is the difference between horizontal and vertical dominoes. These can 
be subdivided further. Colour the underlying lattice squares black and white according to a 
checkerboard fashion in such a way that the left square on the top line is black. Let a horizontal 
domino be of type N or north if its leftmost square is black, and of type S or south otherwise. 
Likewise let a vertical domino be of type W or west if its topmost square is black and type E or 
east otherwise. In figures [T] and [2] the S and E type dominoes have been shaded for convenience. 

One way of sampling from this measure is the so called shuffling algorithm, first described 
in [3], and very nicely explained and generalised in It is an iterative procedure that given 
a random tiling of A n and some number of coin-tosses, produces a random tiling of a diamond 
of A n+ \. You start with the empty tiling on Aq and you repeat this process until you have a 
tiling of the desired size. It is a theorem that this procedure gives all tilings equal probability, 
provided that the coin-tosses we have made along the way are fair. 
The algorithm works in three stages. Start with a tiling A n . 

Destruction: All 2x2 blocks consisting of an S-domino directly above an N-domino are 
removed. Likewise all 2 x 2 blocks of consisting of an E-domino directly to left of a 
W-domino are removed. 
Shuffling: All N, S, E and W-dominoes respectively move one unit length up, down, right 
and left respectively. 



ON THE SHUFFLING ALGORITHM FOR DOMINO TILINGS 



7 



Creation: The result is a tiling of a subset of A n+ i. The empty parts can be covered in a 
unique way by 2x2 squares. Toss a coin to fill these with two horizontal or two vertical 
dominoes with equal probability. 

Figure [T] illustrates the process. In the leftmost column there are tilings of successively larger 
diamonds. From column one to column two, the destruction step is carried out. From there to 
the third column, shuffling is performed. These figures contain several dots which will concern 
us later in this exposition. The creation step of the algorithm applied to a diamond in the last 
column gives (with positive probability) the diamond in the first column on the next row. 

To study more detailed properties of random tilings it is useful to introduce a coordinate 
system suited to the setting and a particle process such that the possible tilings correspond to 
particle configurations. 

In the left picture in figure [H the S and E type dominoes are shaded and a coordinate system 
is imposed on the tiling. For each tile there is exactly one of the x lines and exactly one of the 
y lines that passes through its interior. Indeed we can uniquely specify the location of a tile by 
giving its coordinates (x,y) and type (N, S, E or W). You can see that along the line y = k 
there are exactly k shaded tiles, for y = 1 . . . 8 where 8 is the order of the diamond. The obvious 
generalisation of that statement is true for tilings of A n for any n. We shall call the occurrence 
of a shaded tile a particle. The right picture in figure [21 is the same tiling but with dots marking 
the particles. 

Just to fix some notation, let x\ be the x-coordinate of the i:th particle along the line y = j. 
It is clear from the definitions that these satisfy an interlacing criterion, 



We will now see how the shuffling algorithm described above acts on these particles. 

It turns out that the positions of the particles is uniquely determined before the creation stage 
of the last iteration of the shuffling algorithm, and we have marked these with dots in the last 
column in figure [TJ As can be seen in that figure, running the shuffling algorithm to produce 
tilings of successively larger Aztec diamonds imposes certain dynamics on these particles. That 
is the central object of study in this article. 

Let us first consider the trajectory of x\. As can easily be seen in figure [1] on the y = 1 
line there are always a number of W-dominoes, then the particle, then a number of N-dominoes. 
Depending on whether the creation stage of the algorithm fills the empty space in between these 
with a pair of horizontal or vertical dominoes, either the particle stays or its x-coordinate will 
increase by one in the next step. Thus the first particle performs the simple random walk 



Consider now the particles on row y = 2. For x\ , while x\ (t) < x\ (t) it performs a random walk 
independently of x\, at each time either staying or adding one with equal probability. However, 
when there is equality, x\{t) — x\(t), then the particle must be represented by a vertical (S) 
tile. Thus it does not contribute to growth of the west polar region, thus the particle will remain 
fixed. In order to represent this as a formula, we subtract one if the particle attempts to jump 
past x\. 



(24) 



r 






(26) x\{t) =x\{t-l) + 7l 2 (t) - l{xf(t - 1) + 7l 2 (i) = 
Symmetry completes our analysis of this row with the relation 

(27) x 2 2 {t) = x\(t - 1) + 7 |(t) + \{x\{t - 1) + 7 |(t) 



x\(t-l) + i} 



x\(t-l)}. 




Figure 1. The shuffling procedure. S- and E-type dominoes are shaded. 



ON THE SHUFFLING ALGORITHM FOR DOMINO TILINGS 



9 



For the third row, our previous analysis applies to the first and last particle. 

(28) x\{t) = x\{t - 1) + 7i 3 (i) ~ l{xf(t - 1) + jfit) = x\{t - 1) + 1} 

(29) x\{t) = 4(t - 1) + 7 I(t) + 1{4(* " 1) + 7 3 3 (*) = 4(t - 1)} 

On y = 3 between a;^ and there must be first a sequence of zero or more E dominoes, 
then x%, then a sequence of zero or more N dominoes. While x\ is in the interior of this area it 
performs the customary random walk. It must interact with x\ and x\ in the same way as we 
have seen other particles interacting above. 

So 

4(t) = 4(t - 1) + ii(t) - HxUt - 1) + 7 2 3 W = x 2 2 (t-i) + 1} 

+ l{xl{t-l)+ 1 l{t) = x\{t-l)}. 



The same pattern repeats itself evermore. 



(31) 


4(t) 


= arj(t- 


i)- 






1)- 




7- 

= 




i) + 




(32) 


^(*) 


= x){t- 


i)- 




+ l{a^(t- 


1)- 




7- 


:!(*- 


i)} 




(33) 




= 4(t- 


i)- 




-i{4(t- 


1)- 


H7|(*) 


7- 




i) + 


i} 


(34) 










+ l{xj(t- 


1)- 


WW 


= 4 


:!(*- 


lU- 





with initial conditions x\{j) = i for j = 2, . . . and 1 < i < J. 




Figure 2. Same diamond 



In order to analyse this situation it is suitable to perform a change of variables, 
(35) Xi(t)=x>(t-j) ) 
which gives the equations given in section [2J 

5. Transition probabilities on two lines 

In order to analyse the dynamics just described we follow Warren's example and first consider 
just two lines at a time. What we do in this section is very similar to section 2 of [T2"] . 
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Consider the W n+1,Tl -valued process process (ci n (t)) = (X(t), Y(i)) with components X±(t), . . . , X n 
and Yi(t), . . . , Y n (t), satisfying the equations 

Y i (t + l) = Y i (t)+0 i (t) 

Xi(t + 1) = JCi(t) + ai(t) - l{JTi(t) + ai (t) =Y 1 (t + l) + 1} 

(36) X,{t + 1) = Xi(t) + ai(t) + + on(t) = Yi-^t + 1)} 

-l{X i (t)+a i (t)=Y i (t + l) + l] 
X n+1 (t + 1) = Xi(t) + a n+1 (t) + l{X n+1 (t) + a n+1 (t) = Y n {t + 1)} 

where oti(t) and 0i(t) are i.i.d. coin tosses, s.t. P[a,(t) — 0] = P[a,(t) = 1] = |. They evolve 
until the stopping time r = min{t : Y{(t) — Yi + \(t) for some i S {1, . . . , n — 1}}. At the time r 
the process is killed and remains constant for all time after that. This is a very simple dynamics, 
each Yi either stays or increases one independently of all others. The Xj do the same but are 
sometimes pushed up or down by a Y^\ or Yi respectively so as to stay in the cone ■w n,n+1 . 
This is the discrete analog of the process Q defined in section [3] of this paper. 

Lemma 5.1. For any f : <w n+1 > n — > R, 

(37) J2 JoCC*. y)' y'))/((*', y')) = /((*, v))- 

Proof. Let m = 2n + 1 and z\ = x±, z 2 = y\, . . . , z m -\ — y n , z m = x n +i. Equation (42) in |12j 
states that 

(38) det {^f/^ i< j j\ = ^ <^^<z' 2 ,...,z m <z' m } 

for z, z' £ , W n . Applying the operator A z / (— A a2 )A z / . . . (— A Zm _ 1 )A z ^ to both sides of that 
equality turns the left hand side into io((x,y), (x' ,y')) and the right hand side into l{zi = z{, 
z 2 = Z 2 , . . . , z m — z' m }. □ 

Proposition 5.2. q t , for t = 0, 1, . . . , are the transition probabilities for the process (X, Y), i.e. 
for{x,y),{x',y')€W n+1 ' n , 

(39) fi({x, y), (x',y')) = pC**) [X{t) = x', Y(t) = y';t< r] 
Proof. Take some test function / : •w n+1 ' n — > M. Let 

(40) F(t,(x,y)):= £ tf((x, y), (x', y'))f{x', y>) 
and 

(41) G(t, (x, y)) :=E^[f(X t ,Y t );t<r} 

We want of course to prove that F and G are equal and we will do this by showing that they 
satisfy the same recursion equation with the same boundary values. By lemma 15.11 we already 
know that 

(42) F(0,.) = G(0,.)=/(.). 
The master equation satisfied by G is 

(43) G(t + 1, (x, y)) = ^ ^2 Xl + ai > m + bl ' X2 + ° 2 ' ' ■ ' ' Vn + bn > Xn+1 + a «+i)- 

a 4 ,6ie{0,l} 

This formula simply encodes the dynamics that each particle either stays or jumps forward one 
step. This needs to be supplemented with some boundary conditions that have to do with the 
interactions between particles. 
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When two of the j/i-particles coincide, this corresponds to the event t = t, which does not 
contribute to the expectation in (|TT|) . Thus 

(44) G(t, . . .,yi-i = z,Xi,Vi= z,...) := 0. 
Also, the particle Xi cannot jump past yi, 

(45) G(t, . . .,Xi = z + l,yi = z, .. .) := G(t, . . . ,Xi = z,yt = z, . . .) 
and Xi+i must not drop below g/j + 1, 

(46) G(t, . . . ,yi = z,x i+ i =z,...) := G(t, ...,y l = z,x i+x = z + 1, ...)• 

G(t+1, •) is uniquely determined from G{t 1 •) using the recursion equation and boundary values 
above. It follows that G is uniquely defined by the recursion equation (|43p and the boundary 
conditions (|42I44|45I46I) . 

Observe that, all functions g : Z — > R, satisfy 



(47) 



^(g(x)+g(x + l))=g(x) + ^Ag(x). 



Using this identity many times on (|4T?|) shows that 
(48) 

G(m,(z,y)) = (1+^^(1+^X1+^ 



) Uni 3-n+ 1 ) 



which can be rewritten as 
(49) 



/n+l n \ 

A t G(t,(s,y)) = n( 1 +2 A -)n( 1 +2 A ^)" 1 G (*> (*'»))• 

\i=l i=l / 



The boundary conditions can be rewritten in this notation as well, equations (|44l45l46p can 
be rewritten to 

(50) G(t, (x, y)) = when y { = y i+1 , 

(51) A Xi G(t, (x, y)) — when xi = yi and 

(52) A Xi+1 G(t, (x, y)) = when x l+1 = y 4 . 

Now let us look at F. The observation (|4T|) gives that <f> * i/j = (1 + ^A)?/\ In particular, 

0< t+1 )(xi -xi) A~V (t+1) (yi -a!i) - i </> (t+1) (4 



F(t +l,(a:,»)) 



A^+D^ - yi ) - yi) A^+D(4 - 2/0 

0< t+1 ) (x[ - x 2 ) A" V (t+1) - Xi) (4 - x 2 ) 



' ^(si-si) A-V (t) (yi-a;i)-l </> (t) (4-xi) 
A^+D^-yO 0( t+1 Hyi-2/i) A^ t + 1 )(4- 2/1 ) 
(xi - ar a ) A-V (t+1) (j/i - x 2 ) ^ t+1 ) (x' 2 - x 2 ) 



n+l 



f[(l + lA Xi )f[(l + lA m )F(t,(x,y)) 



i=l i=l 

which shows that F satisfies the same recursion (|4"9")l as G. Now let us take a look at the boundary 
values. 
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qt((x, y), (x' , y')) is zero when yi = yi+\ because two of its rows are then equal. When yi = xi 
for some i then A Xi q"((x, y), (x',y')) — because two rows will be equal when you take the differ- 
ence operator into the determinant. The same argument shows that A. Xi+1 q™((x, y), (x',y')) = 
when yi = Xi+\. Applying this knowledge to the sum F, shows that 

(53) F(t, (x, y)) = when y t = y i+1 

(54) A x .F(t, {x, y)) = when x % = y % 

(55) A Xi+1 F(t, (x, y)) = when x i+ i = y t 

Since F and G satisfy the same recursion equation with the same boundary values, they must 
be equal. □ 

Again, following the example of Warren, we observe that it is possible to condition the pro- 
cesses never to leave •w n > n+1 via a so called Doob /i-transform. See for example [5] for details 
about /i-transforms for discrete processes. Let 

(56) h n (x) = Y[ (xj-Xi). 

l<z<j <n 

The ^-transform of the process above has transition probabilities 

(57) ^ + ((x,y),(x',y')) = (x',y')). 

Just for the sake of notation, call the transformed process (ci n,+ (t)). 

The idea now is to stitch together the process x from processes Q k,+ for fc = 1, . . . , n — 1, 
just like Warren does in the continuous case. For this we need to establish some auxiliary results 
about Q,™ and Q™ ,+ . 

One observation to make is that it is possible to integrate out the x variables from q"' + . 

(58) n> + (y,y') ~ [ £' + {x>y)dx> = detk^ - y^^n 

Jx' 1 <y' 1 <---<y' n <x' n+1 n n\y) 

where dx' is counting measure on •w n+1 = {x G Z™ +1 : x\ < ■ ■ ■ < x n +x}. The reader might 
recognise as a /i-transformed version of the transition probabilities from the Lindstrom- 

Gessel-Viennot theorem. Thus this can be seen as the transition probabilities for a process 
on 'W n , where all n particles perform independent random walks but are conditioned never to 
intersect, i.e. never to leave •W n . We state this as a proposition. 

Fix n > and let x = (1, . . . , n + 1) G Z" +1 and y = (1, . . . , n) G Z n . 

Proposition 5.3. Consider the process (o. n ' + (*)) = (X(t),Y(t)) started in (A(0),F(0)) = 
(x,y). The process (Y(t)) is governed by p n ' + . 

Now for a technical lemma. For x G W™, let •W n (x) — {y G K" : x\ < y\ < ■ ■ ■ < y„ < 
x n+ i} C Z™ and for y G let 

(59) x n {x,y)=n\-^- y 

It is not a difficult calculation to show that X n (x, •) is a probability measure on W n (x). Just 
rewrite h n (y) as a Vandermonde matrix and perform the summation over all y. 

Lemma 5.4. 

(60) f \ n (x, y)q?- + ((x, y), (x' , y')) dy = Pt " +1 ' + (^, x')\ n (x', y') 

JW"(x) 

where dy is counting measure. 
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Proof. An elementary calculation given the explicit formula for q t ' . □ 

Theorem 5.5. Consider the process = (X(t),Y(t)) started in (X(0),F(0)) = (x,y). 

The process (A(i)) is governed by p n+1 ' + . 

Our proof of this is very similar to the proof of proposition 5 in [12j . 
Proof. Recall the the transition probabilities for (o," ,+ (t)) are j,™ . 

¥[x(h) e A u ...,X(t k ) e A k ] = 

dyi--- / dy k q^ + ({x 1 y),(x 1 ,y 1 ))qt^ tl ((x 1 ,y 1 ),(x 2 ,y 2 )) ■ ■ ■ 

■ ■ ■ it*-t h -x(( x k-i>yk-i), {xk,yk)) = 

•W n {x) contains but one element. 



/ dx\ ■ ■ 


■ dx k 




J A k J w 



dxi-- - I dx k I dy dyi ■ ■ ■ dy k 

Ai J A k JW n (x) Jw n (x 1 ) JW n (x k ) 

X n (x,y )q^ + ((x,y ), (x 1 ,y 1 ))q^ tl ((xx,y 1 ), {x 2 ,y 2 )) • ■ ■ 

• ' ' ?t" -i h _! ((xk-i, Vk-i), (xk,Vk)) = 

Repeated applications of lemma 15.41 conclude the proof. 

dyi"- / dy kP ^ + (x,x 1 )X n (x 1 ,y 1 )q^ ti ((x 1 ,y 1 ),(x 2 ,y2)) ■ ■ 

• • • itk-t k ^(( x k-i,yk-i), {xk,yk)) = 

dykp?; + {x,x 1 ) P ^- tl {x 1 ,x 2 ) ■ ■ -Pt k -t k _ 1 (x k -i,x k )X n (x k ,y k ) = 
dxi-- - / dx feF ™; + (5, Xi)p^ h (xi,x 2 ) ■ ■ ■ p£*'-t*_i (x k -i,x k ) 

Ax JA k 



I dx\ ■ ■ 


■ / dx k I 


lAt 


J A k Jw 


I dx\ ■ ■ 


■ dx k 


lAt 


J A k J w 



□ 



6. Transition probabilities for the Aztec Diamond Process 



Let us now return to the process (x(t)) that came from the shuffling algorithm. Observe in 
the recursions ([S]), the formulas that define X k+1 contain X k but not A- 7 for j < k. Thus X k 
is conditionally independent of (A 1 , . . . , X k ~ 1 ) given X k . Also, the dependence of X k+1 on X k 
is the same as the dependence of A on Y in Q_ k and in Q k ' + , see - This, together with 
theorem 15.51 lends itself to an inductive procedure for constructing the process x . 

(1) The process (X k (t),t = 0, 1, . . . ) is started at A fe (0) = x k and has transition probabilities 
governed by p k ' + for k = 1,2, . . . , k, 

(2) By proposition 15.31 X k can be considered as the Y component of the process Q k,+ . By 
the observation above, the pair of processes (X k (t), X k+1 (t)) has the same distribution 
as Q k,+ started at (x k ,x k+1 ) and are thus governed by transition probabilities q k,+ , 

(3) The process A fe+1 is conditionally independent of (A 1 , . . . , X k ~ 1 ) given X k . 

(4) By theorem 15.51 the process X k+1 is governed by transition probabilities p k+1 ' + and 



started at A fe+1 (0) = 
This proves theorem 12.11 



fe+i 
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7. ASYMPTOTICS 

We shall now see some results that lend support to our conjecture that Warren's process X 
can be recovered as a scaling limit of x, the process from the Aztec diamond. Let us rescale time 
by t = Nt and space by Xi = \Nt + ^y/Nxi for i = 0, . . . , n in the above processes. 

First let us show how to recover the entrance law for Warrens process. Recall that the Aztec 
diamond process the discrete process starts at X™(0) = x n and X™ +1 (0) = x n+1 . 

Lemma 7.1. 

2n+l 



(61) ?t "< + ((*",x" +1 ), - v?{x,y) 

as N — > oo where v"(x,y) is given by 119\) . 

Proof. Observe that we can write 
(62) 

2-^VN 2n+ \^ + ((x n+ \x n ),(x,y)) = 2- 2 "+W 2 ' i+1 f X n (x,y) q ^ + ((x,y),(x,y)) dx" 

where dx n is counting measure on the space W (x n+1 ) which has only one element. Then we 
apply lemma [5T31 

(63) = 2- 2n+1 VN 2n+1 p^ + (x n +\x)\ n (x,y) 

which can be written explicitly as 



2"' 



t 

Xi - j 



l<i,j<n+l 



(64) =2 -^VN 2n+1 n\— ^^det 

which can be evaluated by a formula of Krattenthaler (Theorem 26 of [9])- Applying Stirling's 
approximation to the result shows our theorem with Z n = (27r) n / 2 \[. <n j\. □ 

Likewise, Warren's expression for g™ can be recovered as a scaling limit from our expression 
for f™. By Stirling's approximation, 

(65) lVN<f>®(x)^<pt(x), / <pt(y)dy 
and 

(66) l-NA^(x)-, ^ t (x) 

4 dx 

uniformly on compact sets as N — > oo where 

(67) <pt(x) = -^=e~^. 

y Ant 



Lemma 7.2. For Xi = \Ns + hy/Nxi, x\ — |iV"(s + t) + ^yNx'i, and the same relations for y 



and y' , 

(68) {^f ) 2n+ \ n m{{^ V), (x', V')) -> y), (x',y')) 

uniformly on compact sets as N — » oo. 

Proof. Just insert the limit relations given above for <f> and tp in the explicit expression for . □ 
Finally, one of the main results of this article is the following. 
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Theorem 7.3. The process (X k (i),X k+1 (t)) from x , extended by interpolation to non-integer 
times t, converges in the sense of finite dimensional distributions to the process (X k (t), X k+1 (t)) 
from X. 

Proof. For times ti, . . . , t m , and compact sets A±, . . . , A m G \y n > n + 1 ) W e need to study 



(69) lim / dxidiji ■ ■ ■ / dx m dy m x 



x ?£(0,2/)> (S: 1 ,y 1 ))qg_ ii {{x 1 ,y 1 ), (£2,2/2)) • • • ?f m -t m -i ((^ro-i, jfm-i), (x m ,y m )) 

where of course dxidjji is point measure on <w n ' n+ . This is a case of a Riemann-sum converging 
to an integral. By the uniform convergence of in lemmas 17.11 and 17. 2\ we can interchange the 
order of integration and taking the limit. This gives 



(70) = / dxidyi- ■■ dx k dy k x 

J A x J Ak 

x ^"(^l^iX-t! ((xi,yi), (2:2,2/2)) • ■■q? k _ th _ 1 ((x k -i,yk-i), {xk,yk)) 

which proves our theorem. □ 

Theorem 11.11 follows from theorem 17.31 by just restricting to (X k ). We feel that given this 
theorem together with the fact that X k+1 is conditionally independent of X 1 , . . . ,X k ~ 1 given 
X k , lends a lot of credibility to the conjecture in the introduction. 

We can also say something about the limit at a fixed time. Let Ki be the cone of points 
x — (x 1 , . . . , x n ) with x k = [x\, . . . , xfy e Z fe such that 

(71) x k+1 <x k <x k +l 

For each x n € •W n we will denote by x(x n ) the set of all (a; 1 , . . . , a;™ -1 ) such that (x 1 , . . . , x 71-1 , x n ) S 
K_. The number of points in x{x n ) is 

(72) card( ^(^ ))= M^l. 

It follows from the characterisation in section[6]that at a fixed time the distribution of A"™ -1 ^) 
given X n (t) is X n ^ 1 (X n (t), •). Together with the conditional independence noted in that section 
this implies that the distribution of (X 1 ^), . . . ,X n ^ 1 (t)) given X n (t) is uniform in <K_{X n {t)). 
So the probability distribution of X (t) is 

V^T^ T. \ "v(x n ~ T n \ 

(73) m?(x)= P ?(x n ,x n ) X{ ' \'/ / w 
y ' card(^(x")) 

where x( xk > x k+1 ) is one iff x k+1 < x k < x k ^ for alii = 1, . . . , k — 1 and zero otherwise. 

For x n € W n , define K(x n ) as the set of (x 1 , ... , a;™ -1 ) where x k e M fc satisfying a;^ 1 < < 
x i+i- The n(n — l)/2-dimensional volume of K(a?™) is 

,74, »^»^- 
Theorem 7.4. Consider the process (x(t)) t( z%+ under the reseating 

k _ X k {N) - \N 



(75) X 

2 

As N — > 00, X — > A weakly where A /ias distribution 

X{x 1 1 x 2 )... x {x n ~ 1 ,x n ) 



(76) 



vol(K(a; n ) 



16 



ERIC NORDENSTAM 



The expression for /i™ is given in equation ([2"2f . 

Proof. Put in the correct rescaling in (173j) and perform a computation that is practically the 
same as that in the proof of 17.11 □ 

This distribution A also happens to be the distribution of X(l), which is consistent with our 
conjecture. It has been studied pQ and is the distribution of the GUE minor process mentioned 
in the introduction. This proves theorem 1 1.21 

8. Closing Remarks 

Looking at the expression of transition probabilities it is natural to ask the question, 
what happens if we plug in a different <\> than + Si) into that determinantal formula? It 
turns out that for many other <j> this gives a valid transition probability, although we do not fully 
understand why the Doob /i-conditioning still works in that case. It would be interesting to see 
sufficient and necessary conditions on cj) for this construction to work. 

We should mention an article by Dieker and Warren, [2]. They study only the top and 
bottom particles separately from our model, i.e. in our language (Xf (t), Xf (i), . . . , X™(i)) and 
(X{ (t) , X\ (t) , . . . , (t) ) from x (t) . They consider both geometric jumps (<j> = (1 — q) (5q + q8\ + 
q 2 S 2 + ...), for < q < 1) and Bernoulli jumps (</> = pSo + qSi where p + q = 1). They write 
down transition probabilities but do not do the rescaling to obtain a process in continuous time 
and space. 

Another reference worthy of attention is [I], by Johansson. He considers only geometric jumps 
and studies only the top particles (X-J (t) , X| (t) , . . . , X"(i)) from our model, with a slight change 
of variables that is of no real importance. He not only writes down transition probabilities, but 
also recovers the top particles Warren's process X as the limit of his process properly rescaled. 

All the results proved in this article can be generalised to <fi = P$o + q$i where p + q = 1. It 
is also not very difficult given my calculations to write down transition probabilities for the top 
particles and to rescale that to obtain the top particles in Warren's continuous process, analogous 
to Johansson [7] but with Bernoulli jumps. We have not written included that calculation here 
since we don't think it ads much to our knowledge of these processes, but it is a fact that adds 
to the plausability of the conjecture of this article. 
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